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Abstract. Hydrodynamical description of the "Little Bang" in heavy ion 
collisions is surprisingly successful: here we systematically study propagation 
of small perturbations treated hydrodynamically. Using analytic description of 
the expanding fireball known as the "Gubser flow", we proceed to linearized 
equations for perturbations. As all variables are separated and all equations solved 
(semi) analytically, we can collect all the harmonics and reconstruct the complete Green 
function of the problem, even in the viscous case. Applying it to the power spectrum 
we found acoustic minimum at the m = 7 and maximum at m = 9, which remarkably 
have some evidence for both in the data. We estimate effective viscosity and size of the 
perturbation from a fit to power spectrum. The shape of the two-point correlator is 
also reproduced remarcably well. At the end we argue that independent perturbations 
are local, and thus harmonics phases are correlated. 

i 

The actual talk was based on papers [H [21 [31 IH E] and included discussion of the 
"jet edge", correlations between harmonics and other general ideas, in particular the 
relation to Cosmology. Due to lack of space, this written version focuses narrowly on 
the specific results, due to astonisingly good agreement between our calculations and 
wide range of data presented at this historic conference. 

A correct theory usually provides impressive higher order applications; recall e.g. 
Newton's explanation of tides or the higher order QED corrections. The title indicate 
a new round of applications of hydrodynamics, now for small perturbations of the 
exploding fireball. Since this development is new, the first application should be as 
simple as possible. Thus, we only consider (near)central nuclear collisions, in which the 
impact parameter is small enough to be neglected and thus there is no average elliptic 
flow. The zeroth order system (without fluctuations) is thus assumed to be completely 
axially symmetric. We use analytic rather than numerical methods. Therefore, on a 
technical level, we base it on the analytic solution for central collisions of conformal 
plasma objects, to be referred to as the "Gubser flow", see [6]. Obviously here we 
cannot give technical details about this solution: let us just define the comoving 
coordinates p,9 given by: sinhp = — '^^'^ ^^^^ ^ , tan 9 = ^^^2^^f_q2j.2 a combination of 
the transverse distance r and the proper time r = yjt"^ — z'^^ with t, z, being the lab. 
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time and the longitudinal (beam) coordinate. The solution does not depend on the 
longitudinal rapidity and azimuthal angle 0. The parameter q represents the system's 
size. The equation for the linearized first-order perturbations around this solution for 
the temperature perturbation T = Tq{1 + 5) can be written in closed form, for zero 
viscosity it is 

d^5 1 (d^5 1 d5 1 d^5\ 4 . 95 _ 

~ ScoshV [m^ + t^m + ^^d^J ^3% ~ ^ ' 
for longer version with viscosity see These equations allow for separation of variables: 
in fact the bracketed part is nothing else but the angular part of the spherical Laplacian, 
and thus it is solved by the usual spherical harmonics Yim{0,4') and harmonics depend 
on / rather than m. The p part of the equation is treated as the "time dependence". 

Each harmonic can be viewed as an independent oscillator, transferring the energy 
between its potential (pressure) and kinetic (flow) forms, with certain /-dependent 
frequency. We studied the p (time) evolution of the perturbations, with and without 
viscosity and found that the amplitude is greatly reduced for higher harmonics in the 
viscous case, while the oscillation phase is more or less preserved. Both in the Big and 
Little Bangs, the time allocated to the hydrodynamical stage of the evolution is limited 
by the "freezeout time" Tf, after which the collision rates can no longer keep up with 
the expansion. At this time each harmonic has a different phase of its oscillation. In 
the Universe the temperature fluctuation is just read from the sky, and thus nodes of 
Si{fo) correspond to these minima. In the Little Bang one has to calculate the specific 
combination of the temperature and flow perturbations. This includes the calculation 
of how the freezeout surface is modified: but the nodes/maxima of this "observable" 
combination generate the minima/maxima in the power spectrum. This simple physics 
is very robust, the minima/maxima are easily predictable and rather insensitive to many 
details such as dissipation. 

Figjl] shows how this idea works in practice. We assumed that the size of the 
initial pertubation is 0.5 fm, corresponding to the Glauber collision model, and use 
different viscosities as shown below. Our calculation predicts the acoustic minimum at 
m = 7 and the maximum at m = 9. Quite remarkably the very first data on this power 
spectrum, from ATLAS collaboration, show the same minimum and maximum. One 
can further see, that the data points for the largest nonzero harmonics m = 6, 9 are 
incompatible with zero viscosity on the level of many sigmas, and even the ADS/CFT 
value is clearly excluded. The estimated viscosity from this plot is about twice larger 
than the AdS/CFT prediction, namely ^ ^ 2. 

A traditional plot is the two-point correlation function as a function of the relative 
angle. Fig|2] shows its calculated shape, with and without viscosity, to be compared 
with the experimental data. Once again, with an appropriate viscosity, the agreement 
of the shape becomes amazingly good. Note in particular the width of the main peak 
and the hight of a plateau. Moreover, nothing has been fitted to the data: it is the first 
calculation with rather approximate speed of sound and the freezeout time preselected 
in the calculation supeimposed with the first round of data from this conference. (The 
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Figure 1. (Color online) The power spectrum of flow harmonics, versus 
m. The data points are preliminary ATLAS data "S]. Four curves top-to-bottom 
(dashed magenta, dash-dotted red, solid green and dashed blue) are our calculations for 
viscosity-to-entropy ratios Anrj/s — 0, 1, 1.68, 2, respectively. All curves are normalized 
to the m=3 harmonics. 
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Figure 2. The two-pion distribution as a function of azimuthal angle difference A(f>, 
our calculation 4 . Two left plots are for viscosity-to-entropy ratios rj/s — 0, 0.134, 
respectively, while the right one is from ATLAS report [8 . Similar data have been 
reported at this meeting by all RHIC/LHC collaborations. 



delta-function perturbation has been placed at some typical location r = 4.1 fm: how 
it depends on that is discussed in [1] in details.) 

Two opposite scenarios, from incoherent to fully coherent, can be defined as follows: 

(i) The "minimal Gaussian model" is a possibility that the perturbations are in fact 
incoherent noise, fully characterized by the mean squares < 5;(m)^ >. Apparent 
agreement with the Green function calculation appears simply because of weak 
dependence of < 6i{inY > on /. 

(ii) The "maximal coherence models", assumes that harmonics are coherent and the 
physical perturbations are in fact described by our Green function. The physical 
argument for it is that quantum fluctuations in nucleon-nucleon interactions at different 
locations (in the transverse plane) should to a large extent be independent from each 
other. Phase correlations has been demonstrated in [2] based on the Glauber model. 

How can one tell the difference between them experimentally? By going beyond 
the two-point correlators. Indeed, as argued in our paper [2], the three-point 
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correlators allow to measure the "resonance terms" including three harmonics such that 
nil + ^2 + ^3 = 0, and extract the terms with the combinations of phases such as 
cos{miipi + m2ip2 + trisips). For central collisions, for which the two-point correlators 
are ~ 10~^, one would expect the three-point ones to be another factor 30 or so down in 
magnitude. For non-central collisions one may use large value of the ellipticity, provided 
one of the harmonics involved is the second one (e.g. mi = 2) and less statistics be 
needed: but in this case the theory is more complex and yet needs to be developed. 

We end up, reminding the reader, that while the calculations themselves are quite 
technical, the underlying physics is very simple, basically the sound circles similar to 
what one finds while throwing stones into the pond. The sound velocity ~ 1/2 and 
the time till freezeout Tpo ~ 2i? (where R is the nuclear size) , thus the maximal radius 
of the circles (the "sound horizon" ) is simply Hg ~ R. Thus, in terms of the azimuthal 
angle A0 = ±lradian. Strong radial flow dramatically enhances the contrast at the 
edge of the fireball, making small deviations of the freezeout surface into two "horns" 
[Ij , best seen in a specially tuned region of the secondaries, with ~ 2 — 3 GeV where 
the radial flow effects are at its maximum. The circle and horns are well seen in our 
calculation. They were found in hydrodynamical simulations by the Brazilian group 
[7], who went on and provided a simple explanation for the three-peak structure of the 
correlation functions. Indeed, if one calls the two horns + and -, the 4 combinations for 
the particle pairs are + +,- -,+ - and - +. The first two produce a peak near zero relative 
angle and the last two would correspond to the relative angle being twice the (angle 
projected) sound horizon, or about 2 radians. This is indeed what is experimentally 
observed. We thus expect that the relative phases of the harmonics be such as to cancel 
the (non-existing) third horn of the leading m = 3 harmonics. It is still amazing to find, 
that the pictures holds against the data up to the 9-th harmonics! 

Partly an inspiration for this work is clearly the events in Cosmology during the 

last decade, which gave us observations of "the frozen sound" scale, also known as 

the sound horizon, both in the cosmic microwave background (CMB) radiation and in 

the distribution of galaxies. As these observations turn Cosmology into a much more 

quantitative science, we hope their analogues for the "Little Bang" will also fix the 

global parameters in question much better. 
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